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Abstract 

The Feynman checkerboard problem is an interesting path integral approach to 
the Dirac equation in '1+1' dimensions. I compare two approaches reported in the 
literature and show how they may be reconciled. Some physical insights may be gleaned 
from this approach. 
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1 Introduction 

Kauffman and Noyes [5] presented an intriguing derivation of the Feynman checkerboard 
problem[2]. Of particular value, Kauffman and Noyes provide explicit expressions for the 
contributions from the various equivalence classes of paths that contribute to the wave func- 
tion. Explicit evaluation of some simple cases reveals that these expressions are incorrect, 
unfortunately Prior to Kauffman and Noyes, Jacobson and Schulmanjl] sketched a deriva- 
tion of the combinatorial factors that occur in the 'RL' class of paths discussed in Kauffman 
and Noyes [5]. By reverse engineering the arguments given in Jacobson and Schulmanjl], I 
was able to write the correct combinatorial factors for the checkerboard problem as developed 
in Kauffman and Noyes [5]. Section [2] is a precis of Jacobson and Schulman's constructive 
procedurejl]. Section [3] shows how to express these results in the notation of Kauffman and 
Noyes[5j. In Section HJ I show how the propagator for the checkerboard problem may be 
evaluated for the discrete time step case and the continuum limit. The Appendix contains 
a graphical representation of all of the possible paths for the specific case of a 3 x 2 grid so 
that the combinatorial arguments in the main text may be followed constructively. 

2 Jacobson and Schulman's Argument 

For a rectangular grid consisting of r steps of unit length in the 'R' direction and I steps of 
unit length in the 'L' direction, a 'typical' 'RL' path consisting of c corners may be found by 
noting that for an 'RL' path there are exactly 1 + (c — l)/2 left turns and exactly (c — l)/2 
right turns where the first step must be to the right and the last step must be to the left. 
Adding the number of left and right turns together, one sees that l + (c— l)/2 + (c — 1/2) = c 
as it should. Note that, for n steps in a particular direction, there are n + 1 grid points, but 
the first and last grid points are constrained by the equivalence class of the path. Thus, for 
n steps there are n — 1 free grid points at which to place (c — l)/2 turns. As noted above, 
for the 'RL' paths, the last turn must be from right to left so there are (c— l)/2 free turning 
points among r — 1 grid points on the 'R' axis and (c — l)/2 free turning points among I — 1 
grid points on the 'L' axis. These considerations allowed Jacobson and Schulmanjl] to derive 
the following count for the number of 'RL' paths consisting of c corners 



This is related to the expression for Nrr or Nll given by Kauffman and Noyes [5] with some 
important differences. Note the symmetry between I and r in Equation [TJ One may verify by 
explicitly constructing the paths for small dimension grids that the expression for Nrl(c) in 
Equation [1] is in fact the correct one. An example of this procedure for a 3 x 2 grid is given in 
the Appendix. Due to the symmetry in the lower argument of the binomial coefficients, one 
expects that Nrl(c) = Nlr(c). Explicit construction of the allowed paths for small grids 
confirms this. In order to understand this result from a geometric perspective, note that 
each 'LR' path may be inferred from a corresponding 'RL' path by rotating a given 'RL' 
path about the center of the grid by 180 degrees^. The reader can construct a few explicit 

1 Thanks to Kevin Knuth of the UAlbany Physics Department for pointing this out. 
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examples of this procedure by considering the diagrams given in the Appendix. 

The quantity (c — l)/2 is the k index of Kauffman and Noyes[5]. In the notation of 
Kauffman and Noyes [5] 

N RL (k) = N LR (k)=Ct 1) Ct 1) . (2) 

This expression is seen to be closer to the expressions for N RR (k) and Nnik) given in 
Kauffman and Noyes[5J. Equation [2] properly counts the number of allowed paths for both 
the 'RL' and 'LR' classes. It is related to the ipo function of Kauffman and Noyes[5]. Similar 
arguments may be employed to derive the appropriate expressions for N RR (k) = C^^Cj! ^ 

and NiL(k) = Cj? ^C^T^. There is an error in the expressions for N RR (k) and N LL (k) in 
Jacobson and SchulmanJI], although the correction is vanishingly small in the large grid, 
large number of corners limit. The expressions given here may be verified by explicitly 
counting paths for small grids. An example of this procedure is given in the Appendix for a 
3x2 grid. It is seen that N RR (k) is related to the ipR wave function and N LL {k) is related 
to the ip£ wave function of Kauffman and Noyes [5]. 



3 Kauffman and Noyes' Wave Functions 

In order to write a correct wave function which properly accounts for the number of paths 
in a given equivalence class for a grid size r x I, it is only necessary to make the following 
identifications based on the wavef unctions given in Kauffman and Noyes [5] 

*> = £ (-i) (<: - 1,/2 G ( t-;) 2 G ( t;i /2 



oddc>l 



(i-i) 

/2) 



evenc>0 



- 2^ (~ l > C (c/2) °(c/2)-r 



evenc>0 

Here C(c-i)/2 f° r example is a 'generalized binomial coefficient' in the terminology of Kauff- 
man and Noyes [5J. This generalized binomial coefficient may be written as[S] C(c-i)/2 = 

(r-l-A)! 



((c-lV2)!(r-l-A-(c-l)/2)l " 

N 



ote that c ranges over odd values for ip and even values for ip R and ipL- For future 
applications, it is useful to redefine the summation index so that it runs over all non-negative 
integers regardless of class. For ip , substitute (c — l)/2 — > k. For ipi and ip R substitute 
(c/2) — 1 — y k. With this choice of summation index one finds 

^° = ^ 1) \\(r-l~k)\k\(Ll-k)\ (3) 

k>o k\(r-l-k)\(k + !)!(/-! - (k + 1))! U 



xb y r 1)k (r-l)\ 

^ } (jfe + i)!(r-l-(Jfe + l))!jfe!(Z-l-jfe)! 
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As noted above, the summation index k in Equations [3H5] is now over all non-negative 
integers regardless of path class. With the definitions given in Equations [3H5J the following 
derivative identities, related to similar expressions in Kauffman and Noyes [5] may be written 
down. These derivative identities are useful for constructing the discretized version of the 
'1+1' Dirac equation. 



dj> R 
Or 

dipp 
dr 

dipL 
dl 

dl 



00 

~1pR 



In order to compare to the one- dimensional Dirac equation, appropriate linear combinations 
of Equations HTHS are needed. Note that the phase of the Dirac equation given in Kauffman 
and Noyes [5] differs from that used in Jacobson and Schulman[4j. The latter phase convention 
is more commonly used and will also be used here, as this leads to a propagator of the 
following form [I] 

^ a = lim V %(c)(jem ) c (6) 

n. — fc-no ' * 



c>0 



where a, [3 £ {L, R} and c is the number of path changes and n is the number of steps. This is 
the form of the propagator given by Feynman in his original formulation of the checkerboard 
problem^]. This form is chosen for the work presented here in order to facilitate comparison 
with previously published expressions. As Kauffman and Noyes note [5], other choices are 
possible for, e.g., real- valued solutions of the Dirac equation. 
Using the phase convention implied by Equation [61 one finds 



1^2 \ _ ( dipi/dr 
iifit J \ di> 2 /dl 



(7) 



Equation [7] is equivalent to Equation 13 of Kauffman and Noyes [5] with the current phase 
convention. The phase convention implied by Equation [7] leads to the following choices for 
ipi and ip2. 

0i = iipo - ipn (8) 
^2 = #o - i>L- (9) 

Equations Eland IH] are correctly phased with respect to the propagator defined by Equation EJ 



4 Constructing the Propagator 

In order to test whether the phases and wavefunctions chosen here are consistent with pre- 
viously published results, it is useful to perform a consistency check. For this reason, the 
propagator for the checkerboard problem will be derived with the wavefunctions defined 
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by Equations I3HH1 In the limit as the number of corners in the path goes to infinity the 
continuum limit for the propagator may be derived. Observe (cf. Equation [6]) 

K RL = K LR = lim V N RL (c)(iem ) c (10) 

71— >00 ' ' 

c>0 

where c is the number of corners and mo is the particle mass. Note that the path specific 
summation index c is used here instead of the generalized k index. This choice facilitates 
comparison with Jacobson and Schulman's derivation [3]. In a system of units where Ti and 
the speed of light are unity, m has units of 1/length. Recall that for 'RL' or 'LR' paths c is 
odd, N R l gives the number of paths connecting the endpoints of the path with c corners and 
n is the total number of steps in the path. It is useful to define the subsidiary quantity m such 
that r = (n + m)/2 and I = (n — m)/2. Following Jacobson and Schulman's exposition [3], 
note that 

n 2 ( fm\ 2 



rl = T[ 1 -{n))- < u > 

Defing the quantity 7 = a/1 — (m/n) 2 , the product rl in Equation CD] may be written as 
rl = (n/(2 r y)) 2 . The quantity m/n may be interpreted as a dimensionless velocity, as n > 
and m can have either sign. Note that m = corresponds to a particle at rest. If the 
path difference b — a is traversed in a time tb — t a then (m/n) 2 = (b — a) 2 /(tb — t a ) 2 where 
time is measured in light meters, say. Performing time slicing in the usual way, note that 
e — (tb — t a )/n. For a large number of time steps such that c/n the combinatorial 
factors in iV^x,(c) may be evaluated asymptotically. The steps involved are similar to the 
ones used to obtain a Poisson distribution from a binomial distribution [9]. Note that 

N ^ = ( ( c -"l)/2 ) ( (c'-l)/2 ) ' < 12 > 

When r, I 3> c in Equation [T2] the ratios of factorials appearing in the binomial coefficients 
may be approximated as 

(r-l)\ (/-!)! 



N RL (c) 



[(c - l)/2]![r - 1 - (c - l)/2)\ [(c - l)/2]![/ - 1 - (c - l)/2]! 

r (c-l)/2 /(c-l)/2 

[(c-l)/2]![(c-l)/2]! 

(/ r )(c-l)/3 

iiM/f' 



(13) 



The asymptotic propagator using Equation [13] may thus be written 

/ n \ 1 

K RL = (*em ) ^(.emo)- 1 ^- U(r _ um]2 - (14) 



2 7 7 [[(c-l)/2]!]2- 

Using the definition of e and introducing the abbreviation z = mo(tb — t a )/7, the expression 
for the propagator may be put in to the following form 



00 



K RL = (iem Q ) ^(-l) fe (^/2) 2fc /(A:!) 2 (15) 



fc=0 
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where k = (c - l)/2. Note that [I] 

v ; ^ sir (s + a + l) \2J v ; 

where T(s + a + 1) = (s + a)! for integer a such that s + a > 0. Using Equation [161 the 
asymptotic form of the K RL propagator may be written. 

K RL = (iem )J (z) (17) 

A similar procedure can be used to evaluate the K RR and Kll contributions to the propa- 
gator. For K RR paths, one has 

K ** = J™ E ( T c/2 ) ( c/2-1 ) 

where the sum is over all even integers greater than zero. The expression for Kll may be 
found from the substitution r — > I — >■ r. Using the same asymptotic analysis as before, one 
finds 

7 ^ (-l) k (z/2) 2k+1 

where the summation index k = (c/2) — 1. Using the definition of J a in Equation [161 this 
may be rewritten 

1 

K RR « -2rem -Ji(2;). 

n 



Recalling that r = (n + m)/2 and 7 = 1/y 1 — (m/n) 2 , and defining t = ne and x = me, 
where |x| < t, one may show 

t + x 

-»■ -em Ji(z), (18) 

T 

where t 2 = t 2 — x 2 . Making the substitution r — > I — > r one may also write down 

K LL -> -em - — (19) 
r 

Noting that the number of corners is even for 'RR' or 'LL' paths and odd for 'RL' or 'LR' 
paths, one may write down the continuum limit for the propagator[l] by dividing through 
by 2e. Combining Equations [T71 [TS] and [L_H] one may write down the continuum propagator 
in a concise matrix form as follows 

V K LR K LL ) 2 V iJ o(z) —t-Ji{z) ) 

Equation [20] reproduces Jacobson and Schulman's result [3]. This form allows efficient com- 
putation of the two-component wave function at time tt> given the value of the wavefunction 
at time t a . Recall that z = mo(t& — t a )/j. Specifically, 



Mh) 2 I U {z) -^Ji(z) J I M* 



lb) 



(21) 



With this construction, it is clear why the propagator K is decomposed into K RR , K RL , K LR 
and Kll terms. 
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5 Estimating the Number of Significant Paths 



From the discussion given in the appendix, it is clear that the number of paths with a given 
number of corners is not uniform. On a fine grid, it is useful to develop a criterion for 
the paths with a given number of corners that make the most significant contributions to 
the propagator K. The discussion given here parallels that of Jacobson and Schulman [I] 
with the combinatorial factors worked out here. Equation [TH is a useful starting point for 
discussion. Using the definition of z, Equation HH can be rewritten as follows 



K 



RL 



27 



n 



i c (z/2y/[[(c-l)/2}\ 



(22) 



odd c 



We seek that value of c which maximizes (z/2) c /[[(c — l)/2]!] 2 . Following Jacobson and 
Schulman we define exp /(c) = (z/2) c /[[(c - l)/2]!] 2 . Thus 



/(c)=clog(z/2)-21og([( C -l)/2]!). 
Using Stirling's approximation, Equation [231 becomes 

f(c) « clog(z/2) - 2[(c - l)/2] log((c - l)/2) + 2[(c - l)/2] - 



(23) 



(24) 



Setting df/dc — we find Co = 2 where Co is that value of c which extremizes exp/(c). In 
order to verify that Co maximizes /(c), compute d 2 f/dc 2 = —1/cq. In the neighborhood of 
Co, therefore 



exp /(c) « exp /(c ) exp 



(c- c j 
2c 



(25) 



Thus the important contributions to Krl come from a narrow range of c values centered on 
Co = z. Similar calculations may be done for Krr and Kll- The quantity Klr = Krl. The 
results are 



K 



RR 



K RL 

K LR 
K LL 



even c > 
27 



c 2 2 - 



2! 



2/2)! (2/2) 



exp 



odd c> 
Krl 



[z/2)\{z/2) 



exp 



2^ 



C-2 



2^ 



E 



fi) c 2 2 



even c > 



[z/2)\(z/2) 



exp 



2^ 



Note that the summands in this limit differ only by a phase. The work of Ord and col- 
laborators [6l [8] has investigated the implications of the cyclic nature of the propagator 
phase and have used it to suggest a way of extending the Feynman checkerboard problem 
to higher dimensions. It is also interesting to note that there is an alternative method of 
evaluating the propagator, based on the Ising model as discussed by Gersch [3]. Schulman 
[TO] has references to the papers of Ord and coworkers in addition to alternative approaches 
to extending the checkerboard model to higher dimensions. 



7 



Acknowledgements 



KAE thanks the University at Albany for partial support of this work through its Faculty 
Research Award Program. KAE also thanks ACERT (NIH NCRR P41 RR016292) for the 
use of its computational resources. 

References 

M. Abramowitz and I. Stegun. Handbook of Mathematical Functions with Formulas, 
Graphs, and Mathematical Tables. Dover, 1964. 

Richard P. Feynman and A. R. Hibbs. Quantum Mechanics and Integrals. McGraw-Hill, 
1965. 

H. A. Gersch. Feynman's Relativistic Chessboard as an Ising Model. International 
Journal of Theoretical Physics, 20:491-501, 1981. 

Theodore Jacobson and Lawrence S. Schulman. Quantum stochastics: the passage from 
a relativistic to a non-relativistic path integral. J. Phys. A: Math. Gen., 17:375-383, 
1984. 

Louis H. Kauffman and H. Pierre Noyes. Discrete physics and the Dirac equation. 
Physics Letters, A218: 139-146, 1996. 

G. N. Ord. A Reformulation of the Feynman Chessboard Model. Journal of Statistical 
Physics, 66:647-659, 1992. 

G. N. Ord. Classical Analogue of Quantum Phase. International Journal of Theoretical 
Physics, 31:1177-1195, 1992. 

G. N. Ord and D. G. C. McKeon. On the Dirac Equation in 3+1 Dimensions. Annals 
of Physics, 222:244-253, 1993. 

F. Reif. Foundations of Statistical and Thermal Physics. McGraw-Hill, 1965. 
L. S. Schulman. Techniques and Applications of Path Integration. Dover, 2005. 



8 



A Explicit Construction on a 3x2 Grid 

In order to get a feel for the practical aspects of the combinatorial factors used here, consider 
the problem of enumerating all allowed paths with one, two, three or four corners on a 3 x 2 

grid. Note that there are {R+L) C R = ( L J total pathfl With R = 3 and L = 2 there 

are 5!/(3!2!) = 10 possible paths. Here are, respectively, the 'RL' and 'LR' paths with one 



corner 









'RL' 













'LR' 







Here are the 'RR' paths with two corners 









'RR' 













'RR' 







Here is the unique 'LL' path with 2 corners 









'LL' 







Here are the 'RL' paths with 3 corners 









'RL' 













'RL' 







Here are 'LR' paths with three corners 









'LR' 













'LR' 







Here is the unique RR path with 4 corneres 









'RR' 







Note that this exhausts all the allowed paths on a 3 x 2 grid. 



2 I am indebted to Prof. Kevin Knuth of the University at Albany Physics Department for this observation 
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